Quantum Encodings in Spin Systems and Harmonic Oscillators 
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We show that higher-dimensional versions of qubits, or qudits, can be encoded into spin systems 
and into harmonic oscillators, yielding important advantages for quantum computation. Whereas 
qubit-based quantum computation is adequate for analyses of quantum vs classical computation, 
in practice qubits are often realized in higher-dimensional systems by truncating all but two levels, 
thereby reducing the size of the precious Hilbert space. We develop natural qudit gates for universal 
quantum computation, and exploit the entire accessible Hilbert space. Mathematically, we give 
representations of the generalized Pauli group for qudits in coupled spin systems and harmonic 
oscillators, and include analyses of the qubit and the infinite-dimensional limits. 
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Quantum computation may be able to perform certain 
tasks more efficiently than a classical computer; for ex- 
ample, Shor's algorithm Q for factoring prime numbers 
on a quantum computer is exponentially faster than any 
known algorithm on a classical computer. The standard 
model of a quantum computer involves coupling together 
two-level quantum systems (qubits) such that the Hilbert 
space of the system grows exponentially in the number 
of qubits. 

A major obstacle to universal quantum computing 
is the limit on the number of coupled qubits that can 
be achieved in a physical system 0. The use of d- 
dimensional, or qudit, quantum computing enables a 
much more compact and efficient information encoding 
than for qubit computing. Qudit quantum information 
processing employs fewer coupled quantum systems: a 
considerable advantage for the experimental realization 
of quantum computing. The harmonic oscillator is a sys- 
tem that naturally provides qudits as quanta in its en- 
ergy spectrum. Qubits are obtained by restricting the 
dynamics to just two of these quanta, namely the vac- 
uum state |0) and the first excited state |1); e.g., pho- 
tons in cavity QED || and interferometry ph. However, 
the control of entanglement in larger Hilbert spaces is 
now feasible (e.g., orbital angular momentum states of 
photons ||). Our aim is to show that the restriction to 
two-dimensional Hilbert spaces is not necessary and that 
higher-dimensional Hilbert spaces are an advantage, par- 
ticularly when the number of achievable coupled systems 
is limited and entanglement between systems with larger 
Hilbert spaces is physically possible. 

A quantum computer also requires gates, realized as 
the unitary evolution under some Hamiltonian. For 
qubits, a universal set of gates is given by arbitrary SU(2) 
rotations of a single qubit along with some nonlinear cou- 
pling transformation between adjacent qubits generated 
by a two-qubit Hamiltonian . For qudit quantum com- 
putation, the issue of creating a universal set of gates is 



more involved. In particular, it is not possible to treat 
coupled qudits as a collection of qubits, because (typi- 
cally) one does not have access to "pairwise" Hamiltoni- 
ans between two arbitrary levels of coupled qudits. For 
example, in a system of coupled oscillators realized as 
radiation modes in a cavity, Hamiltonians that gener- 
ate single-gate operations such as a coupling of the i th 
level of one oscillator and the j th level of another can- 
not be realized physically. Thus, quantum computation 
with qudits requires an investigation not only into the 
coupling of multilevel systems but also the set of phys- 
ically realizable Hamiltonians with which one can con- 
struct a universal set of gates. In this communication, 
we develop transformations for a collection of coupled d— 
level systems. These transformations are obtained as two 
mathematical realizations of a basis of unitary operators 
for a single qudit. We show how each of these realiza- 
tions can be implemented either in a spin system or a 
harmonic oscillator 0. We establish a SUM gate @, 
which couples qudits and serves as the qudit analogue 
of the controlled NOT gate; this SUM gate employs a 
standard two-mode coupling Hamiltonian. 

The theoretical investigation of qudit computation is 
best expressed in terms of the generalized Pauli group for 
qudits. Recalling the Pauli group for a two-level system, 
a qubit is realized as a state in a two-dimensional Hilbert 
space Ti.2, spanned by two normalized orthogonal states, 
|0) and |1), that serve as a computational basis for Ti.2- 
The unitary operators {X2 = cr x , Zi = er 2 }, where <7i is a 
Pauli spin matrix, generate the Pauli group using matrix 
multiplication: the elements of this group are known as 
Pauli operators and provide a basis of unitary operators 
on Ti.2- 

A qudit is realized as a state in a d-dimensional 
Hilbert space Tid, with a computational basis {|s); s = 
0,1,... , d— 1} serving as the generalization of the binary 
basis {|0), 1 1) } of the qubit. A basis for unitary operators 
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on Tid is given by the generalized Pauli operators [M ^| 
(Xd) a {Z d )\ a,6e0,l,...d-l, (1) 

where A^ and are defined by their action on the com- 
putational basis as follows: 

X d \s) = \s + 1 (mod d)> , (2) 
Z d |s) =exp(27ris/c/)|s). (3) 

The operators Ad and iTd generate the noncommutative 
generalized Pauli group under matrix multiplication, sat- 
isfying 

Z d X d = exp{2ni/ d)X d Z d . (4) 

The analysis of the generalized Pauli group as oper- 
ators in spin systems and harmonic oscillators is neces- 
sary for realizing qudit algorithms and error correcting 
codes 0, [|. For spin systems, we construct the gen- 
erators of the generalized Pauli group in d = 2j + 1 
dimensions using operators that are expressed in terms 
of the SU(2) angular momentum and phase operators. 
This construction allows us to conveniently view a qudit 
as the Hilbert space of a d-dimensional irreducible rep- 
resentation (irrep) of SU(2). For qudits in a harmonic 
oscillator, we obtain a generalized Pauli group generated 
by the number operator N and a phase operator 9. A 
second realization of qudits is given in terms of phase 
states; this realization is "dual" to the first realization 
given here, and allows for the construction of a simple 
SUM gate. By investigating the d — ► oo limit, we show 
that it is not the common generalization of the Pauli 
group for continuous-variable quantum information (i.e., 
the Heisenberg-Weyl group) with position eigenstates as 
the computational basis. 

We begin by constructing a realization of the gen- 
eralized Pauli group for a spin system; i.e., in the d— 
dimensional Hilbert space of a SU(2) irrep of highest 
weight (angular momentum) j = (d— l)/2. Consider the 
standard basis for the su(2) algebra { J z , J± — J x ± iJ y }- 
Let m) z ; m = — j, ... , j} denote the standard weight 
basis for the Hilbert space Hd=2j+i for an SU(2) ir- 
rep of highest weight (angular momentum) j. We use 
a simplifying notation, allowing m to take all the inte- 
ger (or half-integer) values modulo 2 j + 1 , thus defining 

+ 1) 2 = \j, ~j)z- 

With the computational basis defined to be 

\s) = \j,j-s) x , s = 0,l, ... ,d-l, (5) 

the generators of the generalized Pauli group can be ex- 
pressed in terms of operators that act in a natural way 
on SU(2) basis states. Because the basis states are eigen- 
states of J z , we have 

j 

X d ^ \j,m) z (j,m + l\, (6) 
m=—j 

Z d ^eiq>(2iti(j-j z )/d), (7) 



which are unitary and satisfy Eqs. (HH). 

The operators Xd and Zd are conjugate to one another: 

X d = U- 1 -Z d -U, (8) 

where the unitary transformation U is the Fourier trans- 
form in dimension d. It is convenient to view X d as 
the exponent of a Hermitian operator 9 Z , defined so that 
Xd = cxp(27ri 9 z /d), just as Zd is generated by the opera- 
tor J z . The operator 9 Z is known as a phase operator ficfl 
for a spin system. 

The generalized Pauli operators X d and Z& can also 
be realized as operators that act naturally on the space 
Md of dimension d spanned by harmonic oscillator states 
of no more than d — 1 bosons. We define the computa- 
tional basis to be the set of harmonic oscillator energy 
eigenstates 

\s) = \n = s) S o, s = 0,1,... ,d- 1, (9) 

where N\n)no — ^|^)ho- Again, we apply the cyclic 
notation \d) = |0). Now defining the generalized Pauli 
group as operators on this subspace of the harmonic os- 
cillator, the generators Xd and Zd are expressed as 

d-i 

X d Yu I s + < s l > Z d " exp(27riiV/d) , (10) 

s=0 

which are unitary on H^. Again, we view X d as the 
exponent of a Hermitian operator 9 Z , such that Xd = 
exp(27ri 9 z /d); the operator 9 Z is the Pegg-Barnett phase 
operator |ll|, which is well-defined for finite d. We will 
call this representation of the generalized Pauli group the 
number representation. 

An advantage of this explicit realization of X d and Z d 
as unitary operators on the harmonic oscillator Hilbert 
space is that it enables us to explore the d — > oo limit 
in a rigorous way; this limit yields continuous-variable 
quantum computation. The limiting procedure for phase 
operators has been thoroughly investigated |ll], [l^]. In 
this limit, the computational basis remains the harmonic 
oscillator energy eigenstates (now including all states s = 
0, 1, . . . , oo), following Eq. (||). Note that defining the 
phase operator on the infinite-dimensional Hilbert space 
Hoc of the harmonic oscillator presents challenges Jl2| . 
The number operator N and the phase operator 9 Z are 
conjugate in the same sense that momentum and po- 
sition are conjugate, but the limit does not yield the 
usual continuous-variable Pauli group: the Heisenberg- 
Weyl group, with position x and momentum p opera- 
tors as generators. It is also important that the states 
of the computational basis for the limiting case remain 
harmonic oscillator energy eigenstates, rather than posi- 
tion (or momentum) eigenstates or squeezed Gaussians, 
as are commonly used for continuous-variable quantum 
computing. 
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A second realization of X d and Zd in the Hilbert space 
Tid for an irrep of SU(2) can be constructed, with a com- 
putational basis given by SU(2) phase states; this repre- 
sentation is "dual" to the number representation. Con- 
sider the relation 1X2 = exp(i(7r/2)A"2) for qubits; i.e., 
that 



|1) =X 2 \0) = {-iyWQX'lO) 



(11) 



The Pauli operator X2 has two interpretations, each of 
which can be generalized in a different way. In the num- 
ber representation, we interpret X 2 as a cyclic number 
state raising operator |1) — A 2 |0) and generalize this op- 
erator as a cyclic raising operator. However, using the 
relation Jll|), we can also view X 2 as a rotation. (Us- 
ing the su(2) representation X 2 = 2J X , this rotation is 
about the x-axis.) Thus, the state |1) is obtained (up to 
a phase) by rotating |0) by an angle ir about the a;-axis. 
The computational basis states needed for this type of 
generalization to qudits are SU(2) phase states and have 
been investigated by Vourdas (lO) (although using rota- 
tions generated by J z rather than J x ). These states form 
an orthonormal basis for the SU(2) irrep. 

Let {\j, m) x ; m = —j, ... , j} be the weight basis for an 
SU(2) irrep of angular momentum j = (d— 1)/2, where J x 
rather than J z is diagonal; i.e., J x \j, m) x — m\j, m) x . For 
this representation, we define the computational basis 
states to be 



73 EL-j exp(27rims/ d) \j, m) x d odd, 

73 Em=-j exp(27ri(m + \)s/d)\j, m) x d even. 

(12) 



These states form an orthonormal basis for Hd |l(J • They 
are referred to as SU(2) phase states because they are 
eigenstates of a phase operator for spin systems. 

The generalized Pauli operator Xd on this computa- 
tional basis is given by 

X d H-> | eX P( 27ri ^/ rf ) d ° dd ' ( 13 ) 

1 exp(— m/d) exp(27ri J x /d^j d even, 

satisfying Eq. (§). Note that (X d ) d = i for both j in- 
tegral and half integral. The generalized Pauli operator 
Z d is given by 



d-l 



Zd 1— > ^^exp(27ris/d)|s)(s| 



(14) 



which is unitary and satisfies Eq. (|J). Note that we can 
express Zd as the exponent of a Hermitian operator, 



d-l 



Z d = exp(2TTi0 x /d) , § x = ^2s\s){s\; (15) 



s=0 



the operator 9 X is a phase operator for a spin system. 



Note that this representation of the generalized Pauli 
group is "dual" to the number representation of Eqs. (^)- 
(0) in the same sense that the position and momen- 
tum representations of the harmonic oscillator are dual. 
For the number representation, the computational basis 
states are eigenstates of J z , and the phase operator 8 Z 
generates the "ladder" transformations. In the phase rep- 
resentation given here, the computational basis states are 
eigenstates of the phase operator 6 X , i.e., "phase eigen- 
states" , and it is J x which generates the ladder transfor- 
mations via rotations about the x-axis. Both of these 
representations can be considered as natural generaliza- 
tions of the qubit case, because the standard computa- 
tional basis |0) = ||,2)z an( l |1) = l|> — \)z are both 
eigenstates of J z and phase eigenstates of 9 X . 

As with the number representation, this phase repre- 
sentation of the generalized Pauli group can be expressed 
in a harmonic oscillator Hilbert space. Again considering 
the finite Hilbert space H^, the eigenstates of J x (unlike 
J z of the number representation) are replaced with har- 
monic oscillator number states \n) with a boson number 
less than d. The computational basis, then, consists of 
finite-<i phase eigenstates. The generalized Pauli oper- 
ators Xd and Z d are generated by the number operator 
and Pegg-Barnett phase operator, respectively. Again, 
the d — > 00 limit yields challenging problems: it is well 
known that phase eigenstates do not exist in the infinite 
dimensional Hilbert space Hoc of the harmonic oscilla- 
tor 0. 

Despite the issues involving d — > 00 phase operators, 
universal qudit quantum computation is well-defined for 
finite d 0. In the following, we discuss these require- 
ments in terms of an optical realization, where the har- 
monic oscillators are realized as modes in a cavity; such 
a realization has been discussed in 13 . However, this re- 



alization is formally equivalent to any oscillator system. 

To perform arbitrary unitary transformations on a sin- 
gle oscillator efficiently, one may employ a combina- 
tion of linear optics, squeezing, and a nonlinear pro- 
cess such as photon detection or a nonlinear optical 
Kerr interaction [Q. Of particular importance is to re- 
alize the Fourier transform operation on a single qudit, 
which takes number eigenstates to phase eigenstates and 
vice versa. This operation is the generalization of the 
Hadamard transformation for qubits; as it is a unitary 
transformation on a single oscillator, it can be performed 
efficiently as described above. 

For quantum computation, we must also realize a gate 
that performs a two-qudit interaction. A simple con- 
trolled two-qudit interaction gate is the SUM gate 



SUM : \si)i ® |s 2 )2 i-> ki)i ® \si + s 2 (mod d)) 2 



(16) 



Consider two oscillators coupled by the four-wave mix- 
ing interaction Hamiltonian xN\N 2 — X a \ a i a 2 a 2- This 
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Hamiltonian for an optical system describes a four-wave 
mixing process in which \ is proportional to the third- 
order nonlinear susceptibility Let oscillator 1 be 
in a state \s\)\ encoded in the number state basis, and 
let oscillator 2 be in a state s 2 )2 encoded in the phase 
state basis. This interaction Hamiltonian generates the 
transformation 



|*l)x®|(2jf)*i + *2 (mod d)) 2 . 

(17) 



Thus, with fixed interaction time t — 2tt\ 1 , this Hamil- 
tonian generates the SUM transformation on two qudits. 
(Note that a similar gate can be defined for spin systems 
using a J z \ J Z 2 _ type Hamiltonian fllq] .) 

Quantum computation with multiple qudits could be 
performed by coupling several modes in a single cavity; 
each mode realizes a single qudit |u| . Modes are coupled 
via a SUM interaction of the time described above. Note 
that the control qudit for the sum operation must be 
encoded in the number state basis, and the target qudit 
must be in the phase state basis. The encodings of each 
qudit can be swapped (between number and phase state 
bases) using the Fourier transform. 

In summary, we have presented realizations of qudit 
quantum computation in spin systems and harmonic os- 
cillators in terms of number and phase operators. The 
representations of the generalized Pauli group, viewed in 
terms of SU(2) or harmonic oscillator operators, allows 
for qudits to be explicitly encoded into such systems. An 
advantage of this scheme is that the SUM gate employs 
a standard two-mode Hamiltonian to couple two qudits. 
From a rigorous mathematical viewpoint, these realiza- 
tions give natural extensions of the qubit-based Pauli 
group, and allow for the investigation of the d — > oo limit 
and continuous-variable quantum computation. 

By employing qudits rather that qubits, the full size 
of the accessible Hilbert space can be exploited, with the 
advantage of requiring fewer coupled systems for a given 
quantum information process. However, the use of qudits 
requires a different set of quantum gates than the usual 
qubit rotations and two-qubit interactions that are nor- 
mally assumed. The realization of a universal set of gates 
using linear optics, squeezing, and a nonlinear interaction 
is convenient for certain harmonic oscillator systems but 
is not unique; an important challenge is to identify the 



optimal set of gates for a particular system. The analy- 
sis presented here provides the necessary theoretical tools 
for developing qudit quantum computation in spin sys- 
tems and harmonic oscillators as a promising alternative 
to qubit quantum computation. 
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